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The  TOIL  computer  code  is  described  herein  as  it  existed  on 
July  7,  1967.  The  code  has  been  in  continuous  development  for  one  'year 
and  in  its  presented  form  has  been  appljea  successfully  by  General  Atomic 
to  the  kind  of  problems  discussed  later  in  this  report.  However,  the  de¬ 
velopment  and  improvement  of  the  code  are  being  continued,  so  that  duplica-  . 
tion  of  results  (or  even  close  agreement)  between  problems  run  with  the  code 
as  published  and  the  code  as  it  existed  either  before  or  after  this  time 
is  not  necessarily  to  be  expected. 

General  Atonic  has  exercised  due  care  in  preparation,  but  does  not 
warrant  the  merchantability,  accuracy,  and  completeness  of  the  code  or  of 
its  description  contained  herein.  The  complexity  of  this  kind  of  progran 
precludes  any  guarantee  to  that  effect.  Therefore,  any  user  must  make  his 
own  determination  of  the  suitability  of  the  code  for  any  specific  use,  and 
of  the  validity  of  the  information  produced  by  use  of  the  code. 


THIS  DOCUMENT  IS  SUBJECT  TO  SrFCIAl  EXPORT  CONTROLS  AND  EACH 

transmittal  to  foreign  oo>-  '■  o*  fot-tii  imonals  may 

BE  MADE  ONLY  WITH  PRIO"  A  •rr.OiAl  OF  COMMANDING  OFFICER,  U.  S. 
ARMY  ABERDEEN  RESEARCH  «.  DEVELOPMENT  CENTER,  ABERDEEN  EROVINO 
GROUND,  MARYLAND  21005. 


i 


— —  -  A  modification  of  the  OIL  code,  to  treat  two  different  materials  (two 
equations  of  state)  has  been  formulated. 

One  of  the  advantages  that  a  Lagrangian  formulation  has  over  an 
Eulerian  formulation  is  the  ability  and  ease  at  which  material  interfaces 
can  be  represented  and  propagated.  Several  approaches  are  possible  for 
following  a  material  interface  in  an  Eulerian  formulation.  Several  inves¬ 
tigators  are  following  interfaces  (material)  as  they  move  through  the  fixed 
Eulerian  grid.  Others  have  combined  Eulerian  with  Lagrangian  formulation. 

This  scheme  to  be  reported  requires  very  little  modification  to  the 
OIL^  code.  Familiarity  with  the  OIL  code  is  assumed. 

There  are  four  (4)  regions  or  steps  in  the  OIL  code  that  require 
modification  for  treating  two  different  materials; 

(1)  Maintaining  an  interface  between  materials  which  is 
accomplished  by  a  mixed  cell,  in  other  words,  the 
interface  lies  somewhere  in  our  mixed  cell. 

(2)  Calculating  the  pressure  in  a  mixed  cell. 

(3)  The  manner  in  which  the  change  in  specific  internal 
energy  due  to  pressure  forces  only  (PHI)  is  to  be 
distributed  to  a  cell  containing  both  materials. 

(4)  The  scheme  in  PH2  (the  calculation  of  the  transport 
terms)  of  distributing  the  specific  internal  energy  (the 
total  less  the  kinetic)  for  the  cells  containing  both 
materials . 

Item  1  requires  modification  of  the  subroutine  (PH2).  The  discussion 
will  be  concerned  with  the  radial  direction  of  mass  transport,  the  mass  tz’ans 
port  in  the  axial  direction  being  similar.  Calculation  of  the  masses  that 
cross  the  fixed  Eulerian  boundaries  are  performed  as  in  OIL,  however  the 
density  p  in  AM  s  pu  A  At  (u  is  the  weighted  velocity,  A  the  area  and 
At  is  the  time  step)  is  calculated  on  the  basis  of  the  total  mass  in  the 


1 


cell.  Thus,  up  to  this  point,  ue  have  calculated  the  mass  fluxes,  now  we 
must  determine  (for  a  mixed  cell)  how  much  of  each  material  to  move.  Three 
possible  situations  concerned  with  two  materials  arise. 


1.  Material  moving  from  a  non-mixed  cell  to  a  mixed  cell. 
This  presents  no  difficulty  or  modification. 

Example  (non-mixed  to  mixed) 


Mass  flow  is  from  cell  1  to  cell  2. 

AM  =  p1  U  A  At 
AM  =  AM 

X 

AM  .  s  0 

where  =  density  of  x  material  in  cell  1. 

2,  Material  moving  from  a  mixed  cell  to  a  non-mixed  cell  is  cal¬ 
culated  as  follows:  The  acceptor  material  from  the  donor  cell 
is  moved  to  the  acceptor  cell.  If  the  flux  is  such  that  this 
will  more  ■chan  empty  the  acceptor  material  from  the  donor  cell, 
the  excess  is  removed  by  assigning  it  to  the  other  material. 

Example  (mixed  to  non-mixed) 

1  2 

x  . 


Mass  flow  is  from  cell  1  to  cell  2. 

AM  =  p1  U  A  At  where  p1  is  the  total  density  of  both  materials 
in  cell  1  if  AM  >  M^ 

AM  =  M1  ' 

*  • 
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and 


if 


AM 


=  AM 


and 

AMx  =  0  . 

3.  Material  moving  from  a  mixed  cell  to  a  mixed  cell  requires  some 
modification  in  order  to  keep  the  material  interface  defined  in 
a  single  mixed  cell.  The  prescription  or  recipe  if  you  like, 
is  that  each  material  flux  is  weighted  by  the  fraction  of  its 
mass  to  total  in  the  acceptor  cell,  rather  than  the  donor  cell. 

Example  (mixed  to  mixed) 
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Mass  flow  is  from  cell  1  to  cell  2. 

AM  =  p'*'  U  A  At,  where  is  the  total  density  of  both 
materials  in  cell  1,  Then 


and 


AMX 


M2  +  M2 
x 


AM 


M2 

AM  =  -sr-2— p  AM 

•  M  +  it 

x 


Note;  the  superscripts  refer  to  zone  number  and  the  subscript 
to  material  number. 


Item  2  is  the  calculation  of  the  pressure  in  a  mixed  cell.  The 
present  scheme  is  to  iterate  on  the  volume  occupied  by  each  material  until 
the  two  pressures  calculated  are  the  same.  The  partial  volumes  of  each 
material  are  saved  for  use  in  (PHI) , 

Item  3  concerns  the  change  in  specific  internal  energy  for  a  mixed 
cell  due  to  pressure  forces  only  in  (PHI).  The  change  in  internal  energy 
for  the  entire  cell  is 


AQ  =  At  (Vol)  P  ~  • 


The  change  in  specific  internal  energy  for  each  material  is  propor¬ 
tional  to  the  dersity  of  each  material,  or 

AI  -  AS: - - - 

x  f  Mx  M 

T  +  !Uf 


AT  _  _ ^ 

Ai.  “  X^f  M  M 
JS,  + 
f  1-f 

where  (x)  and  ( • )  refer  to  the  two  different  materials  and  f  is  the 
factor  to  multiply  times  the  volume  of  the  total  cell  to  calculate  the  volume 
occupied  by  (x)  material. 

Item  k  deals  with  the  specific  internal  energy  due  to  the  transport 
terms.  The  specific  internal  energy  for  each  material  is  px’oportional  to 
the  specific  total  energy  of  that  material  (the  total  specific  energy  of  the 
cell  plus  that  which  is  transported  in  less  the  amount  transported  out). 


ae  =  m 
x  x 


AE  =  EM 


It.  ♦ 


|  Tfi-y 


Where  the  (~)  tilda  refers  to  the  velocities  and  specific  internal  energy 
after  (PHI). 
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r»  j  1 

Mn  =  A M=  total  mass  of  cell  plus  that  which  is  transported  in  less  the 
x  amo.mt  transpox'ted  out  for  the  x  material. 

m"+1  =  AM  =  similar  term  for  the  dot  material. 

•  • 

Then  the  new  internal  energy  for  the  mixed  cell  is 

«  =  AE.  +  AEx  -  i  ("L.  -  VLl)  +  • 

Thus  each  material  then  has  the  following  specific  internal  energy 

_n+l  *Ex  Q 

lx  =  AM  AE  +  AE 
xx 

Tn+!  _  ^  Q 
.  AM  AE  +  AE 

•  X  • 
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TEST  PROBLEMS 

The  first  test  problem  was  a  comparison  with  SHELL  (a  two-material 
particle-in-cell  code).  A  right  circular  cylindrical  projectile  of  density 
three  times  that  of  the  target  impacted  at  a  velocity  of  1  x  10^  cm/sec 
into  a  semi-infinite  target.  A  separate  equation  of  state  was  used  for 
each  material. 

Figure  1  displays  the  position  of  the  projectile  relative  to  the 
target.  An  (x)  signifies  that  this  cell  was  the  original  projectile,  a 
(dot)  signifies  target  material  and  an  (M)  signifies  that  the  cell  is  mixed 
(contains  both  materials). 

Figure.  2  is  a  similar  plot  for  the  TOIL  results .  The  position  of 
the  material  interface,  as  indicated  by  the  M's,  is  in  good  agreement  be¬ 
tween  the  two  codes. 

The  total  positive  axial  momentum  in  the  grid  as  a  function  of 
time  is  displayed  for  the  two  codes  in  Fig.  3.  The  ~  6$  difference  at 
late  times  is  comparable  to  the  differences  calculated  between  SHELL  and 

oil/2) 

Figures  4  and  5  display  pressure  versus  position  along  the  axis 
into  the  target.  The  position  of  the  shock  front  is  in  good  agreement. 

A  comparison  with  the  one-material  Eulerian  code  (OIL)  was  made 
using  TOIL  for  an  aluminum  right  circular  cylinder  impacting  at  2.6  x  10 
cm/sec  on  a  semi-infinite  aluminum  target.  The  same  density  and  some  equation 
of  state  was  used  for  both  the  projectile  (right  circular  cylinder)  and 
the  semi-infinite  target.  One  would  expect  some  differences  due  to  the 
method  of  transporting  mass  from  a  mixed  to  a  non-mixed  cell. 

Figures  6  and  7  display  the  position  of  the  projectile  relative  to 
the  target.  Figure  8-  displays  the  total  positive  axial  momentum  versus 
time  for  the  comparison. 

Figures  9  and  10  display  pressures  versus  position  along  the  axis 
for  the  two  different  times.  The  agreement  is  very  good,  considering  the 
times  are  slightly  different. 
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A  spherical  hot  source  problem  was  performed  using  two  materials 
with  the  TOIL  cole  and  compared  with  a  spherical  version  of  the  SPUTTER 
(a  one -dimensional  Lagrangian  code). 

Figure  11  indicates  the  original  mixed  interface  between  hot  and 
cold  material,  and  Figs.  12  and  13  display  the  material  interface  as  a 
function "of  time.  The  sphericity  of  the  material  interface  is  indeed  en¬ 
couraging  . 

Figure  14  is  the  pressure  versus  position  (radial  or  axial  from 
TOIL)  for  the  cc.iparison  between  TOIL  and  SPUTTER.  The  near  agreement, 
again,  is  very  encouraging. 

Additional  test  problems  are  being  calculated.  In  addition,  several 
calculations  for  DASA  have  been  successfully  performed  using  the  TOIL  code. 
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0  5  10  15  20 

Fig.  2-~Cell  configuration  foi*  the  TOIL  code  at 
a  time  of  16.-1  jisec 
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Fig.  7— Cell  configuration  for  the  TOIL  code  (Note,  the  distinction 
made  between  the  projectile  and  target  material)  at  t  =  6.19  (Asec 
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Fig.  10 — Pressure  distribution  along 
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Fig.  11.— Cell  configuration  for  the  TOIL  code  (Note,  x  is 
for  the  not  material,  dot  for  the  cold  material  and  M  sig¬ 
nifies  both  materials  on  the  interface)  at  t  =  0 
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Fig.  12--Cell  configuration  for  the  TOIL  code  at  t  =  .96  }Asec 
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Fig.  13  Cell  configuration  for  the  TOIL  code  at  2.25  {1sec 
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